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Abstract. Scale interactions in Hall MHD are studied using both the mean field 
theory derivation of transport coefficients, and direct numerical simulations in three 
space dimensions. In the magnetically dominated regime, the eddy resistivity is 
found to be negative definite, leading to large scale instabilities. A direct cascade of 
the total energy is observed, although as the amplitude of the Hall effect is increased, 
backscatter of magnetic energy to large scales is found, a feature not present in MHD 
flows. The coupling between the magnetic and velocity fields is different than in the 
MHD case, and backscatter of energy from small scale magnetic fields to large scale 
flows is also observed. For the magnetic helicity, a strong quenching of its transfer 
is found. We also discuss non-helical magnetically forced Hall-MHD simulations 
where growth of a large scale magnetic field is observed. 



1. Introduction 

The relevance of two fluid effects has recently been pointed out in several stud- 
ies of astrophysical and laboratory plasmas (Balbus and Tcrqucm, 2001; Sano and 
Stone, 2002; Mirnov et al., 2003; Ding ct al., 2004). The effect of adding the Hall 
current to the dynamics of the flow was studied in several scenarios, particularly dy- 
namo action (Helmis, 1968; Galanti et al., 1995; Mininni et al., 2002, 2003a, 2005b) 
and reconnection (Birn et al., 2001; Shay et al., 2001; Wang et al., 2001; Morales 
et al., 2005). Several of these works showed that the Hall currents increase the re- 
connection rate of magnetic field lines. However, most of the studies of magnetic 
reconnection were done for particular configurations of current sheets. It was shown 
in particular by Smith ct al. (2004) that when a turbulent background is present the 
reconnection rate is dominated by the amplitude of the turbulent fluctuations. The 
process of magnetic reconnection is relevant in several astrophysical and geophys- 
ical scenarios, such as the magnetopause, the magnctotail, the solar atmosphere, 
or the interplanetary and interstellar medium. Reconnection can also play a role 
in the generation of large scale magnetic fields by dynamo action Zeldovich et al. 
(1983). 

Some of the works in Hall-magnetohydrodynamics (Hall-MHD) present conflict- 
ing results, indicating in some cases that the Hall effect can help the growth of a 
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large scale magnetic field (Mininni et al., 2005b) or a large scale self-organization 
process (Mahajan and Yoshida, 1998; Numata et al., 2004; Ohsaki, 2005), while in 
other cases the Hall currents were observed to generate small scales and filamcnta- 
tion (Laveder et al., 2002a, b; Rhcinhardt and Geppert, 2002). 

As a result, it becomes of interest to study the physical processes leading to cas- 
cades and transfer of ideal invariants in three-dimensional Hall-MHD turbulence. 
Phenomena observed in the laboratory and space plasmas tend to show an intermit- 
tent or impulsive behavior (Bhattacharjee et al., 1999) characteristic of turbulent 
flows. The relevance of Hall-MHD turbulence in the solar wind was shown by Ghosh 
et al. (1996). Also, Hall-MHD turbulence can play a crucial role in the transfer of 
matter in the magnetopause as was pointed by Rezeau and Belmont (2001). 

In this work, we study both analytically and numerically three dimensional Hall- 
MHD turbulence as the result of a dynamo process, and from a purely electromotive 
forcing. Detailed studies of shell-to-shell energy transfer from direct numerical sim- 
ulations (DNS) have been done for hydrodynamic (Domaradzki and Rogallo, 1990; 
Ohkitani and Kida, 1992; Zhou, 1993; Yeung et al., 1995; Alexakis et al., 2005a) and 
magnctohydrodynamic flows (Dcbliquy et al., 2005; Alexakis et al., 2005b; Mininni 
et al., 2005a). To the best of our knowledge, the energy transfer in Hall-MHD 
turbulence has not been studied before. 

We show evidence of non-locality of the transfer in Fourier space, and that the 
Hall effect can increase both the transfer of magnetic energy to smaller scales (lo- 
cally), as well as give a novel non-local backscatter of magnetic energy to large 
scales. These results become clear when examining the modification to the tur- 
bulent magnetic diffusivity due to the Hall term. Also, we observe that the Hall 
currents impact on the coupling between the magnetic and velocity fields. The 
transfer of energy between these two fields is different than in the MHD case. The 
Hall-MHD equations also display a backscatter of energy from small scale magnetic 
fluctuations to the large scale flows themselves. The transfer of helicity is briefly 
discussed as well and observed to be quenched by the Hall effect. 

The structure of the paper is as follows. In Sec. 2 we introduce the Hall-MHD 
equations and we define the various transfer terms. In Sec. 3 we derive turbulent 
transport coefficients for the Hall-MHD induction equation. In Sec. 4 we briefly 
discuss the code and details of the mechanically forced numerical simulations for 
completeness. Section 5 presents the transfer terms in Hall-MHD as obtained from 
the numerical simulations. Section 6 shows backscatter of magnetic energy in non- 
helical Hall-MHD magnetically forced simulations. Finally, Sec. 7 summarizes the 
results and discusses implications of our work for the understanding of turbulence, 
dynamo action, and reconncction in Hall-MHD. 



2. The Hall-MHD equations and transfer terms 

In dimensionless Alfvenic units, the Hall-MHD equations are 

dtU + XJ- VU= -VV + B- VB + ^V 2 U + f, (2.1) 

d t B = V x [(U - eJ) x B] + ??V 2 B, (2.2) 

where U is the bulk velocity field, B is the magnetic field, J = V x B is the 
current density, V is the pressure, v is the kinematic viscosity, and r\ is the magnetic 



Energy transfer in Hall-MHD turbulence 



3 



diffusivity. From the Maxwell equations and incompressibility of the flow, V • U = 
V-B = 0. 

The Hall term eJxB in Eq. (2.2) measures the velocity difference between species, 
where the electron velocity is U e = U eJ. Here, e measures the relative strength 
of the Hall effect, with the Hall term being dominant for wavcnumbers larger than 
^Haii ~ 1/e if equipartition between the fields is assumed. The measure of strength 
of the Hall effect can be written as e = Lnaii /Lq where Lq is a characteristic length 
(we will use Lq = 2tt, the size of the box in our simulations). In terms of physical 
parameters, and for a fully ionized plasma, the Hall length is LHaii = cUa/ (ojpiUo), 
where Ua is the Alfvcnic speed, Uq is a characteristic speed, c is the speed of light, 
and LVpi is the ion plasma frequency (when Uq = Ua, Lho.ii reduces to the ion skin 
depth). In a partially ionized plasma, expressions for Lho.ii can be found in Sano 
and Stone (2002) and Mininni et al. (2003a). 

Of special interest is the ratio between the integral length L, the Hall length 
Lho.ii, and the Ohmic dissipation length L v . For LhoII <C L v (e — > 0), the Hall-MHD 
equations reduce to the well known MHD case. In several astrophysical problems, 
such as accretion disks, protoplanetary disks, or the magnetopause (see e.g. Birn 
et al., 2001; Balbus and Terquem, 2001; Sano and Stone, 2002), the Hall scale is 
larger than Ohmic scales although smaller than the integral scale L of the flow. 
We will be interested in this regime in this work, although we remark that the 
separation between these scales in astrophysical or geophysical problems is far from 
what can be achieved in numerical simulations. 

The Hall-MHD equations have three ideal invariants (Turner, 1986). In this work 
we will focus on two invariants, the total energy 



and the magnetic helicity 



where A is the vector potential, V x A = B. These quantities are also ideal invari- 
ants of the MHD equations (e = 0). The third MHD invariant, the cross helicity, 
is replaced in Hall-MHD by the hybrid helicity (Turner, 1986) and is small in the 
simulations we will discuss. 

The expressions we will use for the shell-to-shell energy transfers have been de- 
rived for the MHD case by Verma (2004); Debliquy et al. (2005); and Alexakis 
et al. (2005b). Here we present the derivation of the transfer terms for the Hall- 
MHD equations. Equation (2.2) can be rewritten as 

d t B + U • VB = B • VU - eV x (J x B) + ?7V 2 B. (2.5) 

We introduce a filter in shells in Fourier space, such as Fj< which denotes the 
components of the field with wavcnumbers between K and K + 1 [i.e. Fk (x) = 
Y^k=K F(k)e lk ' x ], from Eqs. (2.1) and (2.5) we can write detailed balance equations 
for the energy, 



d t E v {K) = iy\ -U* -(U • V) ■ U + U A -(B ■ V) ■ B £ 
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- vV 2 XJ K +f-U K \ dx 3 , 



(2.6) 



d t E B {K) = / |X! [- B *"( U ' V) • B Q + B K (B • V) • U Q + 

+ eJ K -(B x J Q )] -77V 2 B k | dx 3 . (2.7) 

Here, Ejj{K) and Eb{K) denote respectively the kinetic and magnetic energy in 
the shell K . The above equations can be written in the more compact form: 

d t Eu(K) = Y^[Tuu(K, Q) + T BU (K, Q)] - vV v (K) + T(K), (2.8) 
Q 

d t E B (K) = Y}TuB(K,Q)+r BB {K,Q))-nV B {K). (2.9) 
Q 

The functions Tuu(K,Q), Tu B (K 7 Q), T BB (K,Q), and T B u{K,Q) express the en- 
ergy transfer between different fields and shells, 

T VU (K,Q) = - y"u K (U- V)U Q dx 3 , (2.10) 

r UB (K,Q) = J\J K (B- V)B q dx 3 , (2.11) 

T BU (K,Q) = Jb k (B- V)U q dx 3 . (2.12) 

In general, for positive transfer, the first subindex denotes the field that receives 
energy, the second subindex the field that gives energy. The first wavenumber corre- 
sponds to the field receiving energy, and the second wavenumber to the field giving 
energy As an example, positive Tjjjj(K,Q) represents energy transfcred from the 
velocity field at the shell K to velocity field at the shell Q. In the same way, positive 
Tjj B {K, Q) represents energy transfcred from the magnetic field at wavenumbers K 
to the velocity field at wavenumbers Q. 

The transfer of magnetic to magnetic energy T BB (K,Q) in Hall-MHD consists 
of two terms 

T BB (K, Q) = TH HD (K, Q) + T B H #V, Q) (2.13) 

TH HI V, Q) = - j B K (U • V)B Q dx 3 , (2.14) 

is the usual MHD transfer of magnetic energy through advection by the bulk velocity 
field, and 

T^ B al \K, Q) = eJj K -(Bx J Q ) dx 3 , (2.15) 

is the transfer of magnetic energy due to the Hall current. Note that the definition of 
the transfer terms corresponds to the MHD case in Alexakis et al. (2005b), except 
for the new term T^^^^X, Q). However, as will be shown later, the behavior of 
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the rest of the transfer terms in Hall-MHD will also be indirectly modified by the 
presence of the Hall effect. 

All these transfer functions satisfy the identity 

T VW (K,Q) = -T WV (K,Q), (2.16) 

where v, w can be cither U or B. This detailed conservation is what allows us to 
define the terms as transfers of energy between shells. Note that other groupings 
of the nonlinear terms in the Hall-MHD equations would not satisfy this symmetry 
condition. 

In Eqs. (2.8-2.9) we also have two dissipation functions and the energy injection 
rate 



vV v (K) ee v J |VU^| 2 dx 3 , (2.17) 
V V B (K) ee nj |VB^| 2 dx 3 , (2.18) 

T{K) ee J f -\J K dx?. (2.19) 

Finally, we can also define the transfer of magnetic helicity. From Eq. (2.2) we 
have 

d t H(K) = ]T T H (K, Q) - rfD H (K), (2.20) 
Q 

where the transfer of magnetic helicity from the wavenumber K to the wavenumbcr 
Q is given by 

T H (K, Q) ee J B K ■ [(U - eJ) x B Q ] dx 3 . (2.21) 

This transfer function satisfies the relation Th(K, Q) = —Th(K, Q). The first term 
in Eq. (2.21) proportional to U is the usual transfer of Hm in MHD, while the 
second term proportional to eJ is the contribution due to the Hall effect. Note that 
as a whole, magnetic helicity is transfercd between the shells K and Q interacting 
with the electron velocity field U — eJ. This is in agreement with the fact that in 
the ideal limit the magnetic field in the Hall-MHD system is frozen to the electron 
velocity field, instead of the bulk velocity field of the plasma as in MHD. 
The dissipation rate of magnetic helicity at the wavenumber K is given by 

rp H {K) = V Jb k - J K dx 3 . (2.22) 

It is also worth noting that, since the magnetic helicity is not a positive defined 
quantity contrary to the energy, the interpretation of its transfer is more difficult. 
We will not attempt here a separation of its different sign components (see e.g. 
Waleffe, 1991; Chen et al., 2003a, b) for the case of kinetic helicity in hydrodynamic 
turbulence) . 



3. Transport coefficients 

In Mininni ct al. (2002), the expression of the a dynamo coefficient was derived 
for Hall-MHD. Although an expression of the Hall-MHD turbulent diffusivity was 
derived by Mininni et al. (2003b), the closure was only valid for specific solutions 
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of the Hall-MHD equations. To interpret the results from the energy transfer, it 
will be useful to have expressions for all the turbulent transport coefficients in the 
induction equation. To this end, and for the sake of simplicity, we will use mean 
field theory (MFT) (Stcenbeck ct al., 1966; Krausc and Racdler, 1980) and the 
reduced smooth approximation (RSA) (Blackman and Field, 1999). RSA was in- 
troduced to solve some ambiguities present in MFT when the magnetic field is 
strong enough to affect the velocity field trough the Lorentz force. Although there 
are still assumptions in MFT not completely justified, at least a qualitative agree- 
ment has been observed with simulations in the MHD case (Brandenburg, 2001) 
and the Hall-MHD case (Mininni ct al., 2003a, 2005b). The transport coefficients 
can also be derived using more elaborate closures, such as the Lagrangian History 
Direct Interaction Approximation (LHDIA) or the Eddy Damped Quasi Normal 
Markovian (EDQNM) closures (see e.g. Lesieur, 1997). It is worth noting that the 
analysis that follows in Sections 4 and 5 is of general validity and independent of 
the assumptions we will use here to derive the turbulent transport coefficients. 
We split the fields into 

U = U + u + u , (3.1) 
B = B + b + b , (3.2) 

where Uo and bo are isotropic and homogeneous solutions of Eqs. (2.1) and (2.5) 
in the absence of the mean fields U and B. The fields with overbars are large scale 
fields, and u and b are small scale corrections to the isotropic and homogeneous 
solutions due to the presence of the large scale fields. The fluctuating fields satisfy 
(u) = (u ) = (b) = (b ) = 0, where the brackets denote an average that satisfies 
Taylor's hypothesis (Krause and Racdler, 1980). Replacing in Eq. (2.5), using the 
equations for the Uo and bo fields, dropping terms quadratic in the fluctuating fields 
u and b, and averaging leads to 

d t B = V x [(U-eJ) xB + e] + r?V 2 B, (3.3) 

where e is the mean field electromotive force 

e = (u e x b + u e x b ) . (3.4) 

Our main aim in this section is to close Eq. (3.3) and write e only as a function 
of averages of the fields Uo, bo, and spatial derivatives of B. A simple argument of 
symmetry shows that in the approximately isotropic case 

s = aB-(3V x B + 7 V x V x B. (3.5) 

From Eqs. (2.1) and (2.5), and subtracting the equations for the mean flows, we 
can also write equations for the evolution of the turbulent fluctuations u and b. We 
drop terms quadratic in u and b, and keep only terms to zeroth and linear order 
in B, 

d t b = V x (U x b + u e x B + eb^ x J + u e x b + u L x b - e) + ?7V 2 b. (3.6) 

In this equation, e involves averaged quantities, and from the Taylor's hypothesis 
it gives no contribution to the mean electromotive force. The fourth and fifth terms 
on the r.h.s. can be dropped using RSA, namely that |u|, |b| <C |B| (note that this 
condition is less stringent than the usual assumptions in MFT, since the amplitude 
of the fields uo, bo can be much larger than the amplitude of the mean magnetic 
field). We will assume the viscosity and diffusivity are small, and as a result we will 
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also drop the last term on the right hand side of Eq. (3.6). Terms proportional to 
U can be removed in the proper frame of reference. Finally we obtain 

dtb pa V x (ug x B + eb e Q x J) . (3.7) 

The second term on the right hand side of Eq. (3.7) involves only spatial derivatives 
of B and gives no contribution to the a coefficient, but is retained here since it will 
give contributions to (3 and 7. 

Following the same steps, we can also write an equation for the evolution of u, 

d t u pa B • Vb + b • VB - Vp. (3.8) 

To obtain the mean field electromotive force we replace time derivatives in Eqs. 
(3.6) and (3.8) by the inverse of a correlation time r. This step, common in MFT, 
assumes the existence of a finite correlation time. At present there is no evidence of 
its validity in general (see e.g. Gruzinov and Diamond, 1995; Blackman and Field, 
2002; Brandenburg and Subramanian, 2005). Since we are introducing a correlation 
time to close these equations, the expressions obtained for the turbulent transport 
coefficients will be considered as symbolic expressions. 

Before replacing the expression for u in Eq. (3.5), Eq. (3.8) has to be solved 
for the small scale pressure p. We will use a technique developed by Gruzinov 
and Diamond (1995) (see also Blackman and Field, 2002). The a effect is linear 
in B and therefore the correct result can be obtained assuming B uniform. Then 
d t u pa B • Vbo. Replacing the time derivative by r^ 1 and replacing the expression 
in Eq. (3.5), we obtain in the weak isotropic case 

a = \ (-Uq ■ V x u e + b • V x b - eb • V x V x ug) . (3.9) 

To compute (3 and 7 we have to keep spatial derivatives of B in Eq. (3.8), and 
therefore we have to solve for the pressure. This was done by Gruzinov and Diamond 
(1995) transforming Eq. (3.8) to Fourier space, and doing a Taylor expansion of the 
projector operator for incompressible u assuming a large scale separation between 
the mean and fluctuating fields. In three spatial dimensions, it was shown that the 
pressure and bo- VB terms give no contribution to e in Eq. (3.4). As a result, we are 
only left with the terms proportional to spatial derivatives of B when Eq. (3.7) and 
the first term on the r.h.s. of Eq. (3.8) are replaced on Eq. (3.4). Again, assuming 
weak isotropy, we obtain the expressions for the remaining turbulent transport 
coefficients 

0=1 (uf + e (uo • V x b^ + b • V x UQ 1 )) , (3.10) 

7 = -y (boX). (3.11) 

The two last terms in a, and the third term in (3 come from the small scale 
momentum equation and are related with the backreaction of the magnetic field into 
the velocity field. In the kinematic regime of a dynamo, a = — r/3 (uq • V x ug), 
which for e = reduces to the MHD case (Krause and Raedler, 1980). The general 
expression for e = reduces to the MHD expression first found using the EDQNM 
closure by Pouquet et al. (1976), a = r/3 (— Uo • V x Uo + b • V x b ). Note also 
that in the MHD case in three dimensions, the turbulent diffusivity (3 = — r (uo 2 ) /3 
is not changed during the nonlinear saturation (Gruzinov and Diamond, 1995). 
The turbulent diffusivity (3 in Hall-MHD is not positive definite, contrary to the 
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pure MHD case (note that negative effective diffusivities can be found in MHD if the 
assumption of homogeneity is dropped, see e.g. Lanottc et al., f999). A negative 
value of (3 represents non-local transfer of energy from the small scale turbulent 
fields to the large scale magnetic field. This result will be of interest in the following 
sections. 

It is worth studying the values of (3 for particular cases. If e is large enough and 
the system is magnetically dominated (Eb 3> Ejj), then (3 f=a —re 2 (j 2 /) /3, where 
jo = V x b and we assumed the average is a spatial average. In this case, (3 is 
always negative implying transfer of energy from the small scales to the large. 

The normal modes of the Hall-MHD equations are circular polarized (V x ug = 
i/sug, V x bo = ±fcbo) and dispersive, and in the limit k ^> 1 they satisfy dis- 
persion relations ui ~ ek 2 B (whistlers, right-handed polarized) and u> ~ B/e (ion- 
cyclotron waves, left-handed polarized). Also, for these waves the fields are related 
by bo = — kB/ui\iQ (Mahajan et al., 2005a). If we assume a background of waves, 
(3 ~ t/3 (uq 2 ) (1 ± 2ek 2 B/u>), which for k large enough can give positive or negative 
turbulent diffusivity according to the orientation of the wave. Note that from the 
dispersion relations, at small scales whistlers give a finite contribution to the tur- 
bulent diffusivity, while ion-cyclotron waves give a much larger turbulent diffusivity 
that grows as k 2 . 

4. Simulations 

In this section we summarize the simulations that will be used to compute the 
energy and hclicity transfer functions defined in Sec. 2. We performed three simu- 
lations in three dimensions with periodic boundary conditions, using a pseudospec- 
tral Hall-MHD code as described in Mininni et al. (2003a, 2005b). Runge-Kutta 
of second order is used to evolve the system of Eqs. (2.1) and (2.2). To ensure the 
divergence- free condition for the magnetic field, a curl is removed from Eq. (2.2) and 
the equation for the vector potential is instead solved, with the Couloumb's gauge 
V • A = 0. The three simulations are done with a spatial resolution of N 3 = 256 3 
grid points. The 2/3 dealiasing rule is used, and as a result the maximum wavenum- 
ber resolved by the code is k m ax = N/3 ~ 85. The kinematic viscosity and magnetic 
diffusivity are set to v = rj = 2 x 10~ 3 , and all the simulations are well resolved, in 
the sense that the kinetic [k v = ((uj 2 ) /v 2 ) 1 ^} and magnetic [k rj = ((J 2 ) /if) 1 / 4 } 
dissipation wavenumbcrs arc smaller than k max at all times. 

In Hall-MHD, the Courant-Fricdrichs-Levy (CFL) condition is more stringent 
than for MHD for which, with equipartition of kinetic and magnetic energy, the 
CFL condition for explicit time-stepping imposes an upper boundary on the time 
step At < Ax/Ua where Ax is the spatial step. In Hall-MHD, the dispersive nature 
of the whistlers impose At < Ax 2 /(eU a)- As a result, smaller time steps will be 
needed as e is increased. Also, since the time step decreases quadratically as the 
spatial resolution is linearly increased, we cannot achieve now spatial resolutions 
higher than 256 3 because of these constraints. 

A helical forcing at fco = 2 given by an ABC flow 

f = [C sin(koz) + B cos(k y)] x + [A sin(fc ^) + C cos(fco^)] V + 

+ [Bsm(k y) + Acos(k x)] z, (4.1) 

with A = 0.9, B = 1, and C = 1.1 was applied in the momentum equation. This 
election of the amplitude coefficients was done to ensure breaking the symmetries 
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Figure 1. Kinetic (thick curves) and magnetic energy (thin curves) as a function of 
time, for runs with e = (solid lines), e = 0.05 (dotted lines), and e = 0.1 (dashed lines). 

of the ABC flow and ensuring a faster development of turbulence (Archontis et al., 
2003). After a first hydrodynamic run made to reach a turbulent steady state, a 
random and small magnetic field was introduced at small scales. Initially the ratio 
of kinetic to magnetic energy was E u /E\, ~ 10 -3 . 

The simulation was continued to see exponential growth of the magnetic energy 
(in the following, we will refer to this stage as the kinematic regime), and finally 
nonlinear saturation of the small scale magnetic field (in the following, Hall-MHD 
turbulence). Three simulations were done, with e = (MHD), e = 0.05 (which 
corresponds to k Hali « 20), and e = 0.1 {k HaU rts 10). Figure 1 shows the time 
history of the kinetic and magnetic energies for these three runs. 

After t w 20, the small scale magnetic fields have reached saturation for all values 
of e, while the large scale magnetic field keeps growing slowly. As e is increased, 
the magnetic energy reached by the system after the non-linear saturation of the 
small scales increases. However, this behavior is not monotonical in e as shown by 
Mininni et al. (2003a). 

The saturation of the large scale magnetic field takes place in a longer time 
(Brandenburg, 2001; Mininni et al., 2005b). Note that one of the biggest challenges 
for DNS is to attain scale separation between the different dynamical ranges that 
must be resolved. Reynolds numbers in simulations are much smaller than the values 
observed in astrophysics and geophysics. Moreover, compared with hydrodynamics 
and MHD, the extra characteristic length scale in Hall-MHD (the Hall scale) makes 
it even harder to achieve a proper separation between all these scales. As a result, 
we will focus in this work in the energy transfer at scales smaller than k (the 
energy injection band), and the late time large-scale evolution of these runs will 
not be discussed here (more details can be found e.g. in Mininni et al., 2005b). 

Figure 2 shows the time evolution of the magnetic and kinetic energy spectra, for 
the runs with e = 0.05 and e — 0.1. As previously mentioned (t « 40) the spectrum 
of energy at scales smaller than ko has saturated and reached a steady state, while 
the magnetic energy at k = 1 keeps growing slowly. Note that the ratio of kinetic 
to magnetic energy at small scales in the saturated state depends on the value of e. 
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Figure 2. Spectra of kinetic energy (thick curves) and magnetic energy (thin curves) as 
a function of time, for t = 5 (solid), 15 (dash-dotted), 30 (dotted), and 45 (dashed), (a) 
e = 0.05, and (b) e = 0.1. 



Another quantity that will be of interest in the next section is the magnetic 
helicity. Figure 3 shows the time history of the absolute value of magnetic helicity 
for the three runs. Note that while in the MHD run (e = 0) the magnetic helicity 
grows monotonically with time, in the Hall-MHD runs the time evolution is strongly 
modified. For e = 0.05 the magnetic helicity grows slower than in the MHD case, 
and for e = 0.1 it changes sign at t ss 24. As was observed by Mininni et al. (2003a), 
the Hall effect inhibits the generation of net magnetic helicity at large scales by the 
helical dynamo process. This inhibition grows monotonically with the amplitude of 
Hall term, and for values of e large enough the magnetic helicity fluctuates around 
zero. The reason for this behavior will be discussed in the next section. 
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Figure 3. Absolute value of the magnetic helicity as a function of time. Labels are as in 
Fig. 1. For e = 0.1 the magnetic helicity changes sign from negative to positive at t ~ 24. 
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Figure 4. Transfer of kinetic energy from Q = 20 to K, Tuu{K, Q = 20), in the kinematic 
regime (solid line), at t = 26 (dotted line), and at t = 45 (dashed line) in the run with 
e = 0.1. 

5. Transfers 

In this section we discuss the energy transfer terms denned in Sec. 2 as obtained 
from the three DNS discussed in the previous section. 

5.1. The run with e = 0.1 

We start discussing in detail the transfer in the Hall-MHD run with e = 0.1. At late 
times in this simulation, when the system is close to cquipartition (Ua ~ U), the 
Hall wavenumber is knaii ~ 10. Since we consider transfer functions at different 
times, for the sake of comparison and unless explicitly said, all transfers in this 
subsection will be normalized using the r.m.s. velocity and magnetic field according 
to their expressions [Eqs. (2.10-2.15)]. Note that since eJ has units of velocity (and 
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U — eJ is the electron velocity), the transfer function T^j^ is normalized using 
^£? 2 |U|). This election also allows for a direct comparison of this term against 

T|| HD [see Eq. (2.14)]. 

Figure 4 shows the transfer of kinetic energy from the shell Q = 20 to kinetic en- 
ergy in shells K for three different times. As previously mentioned, positive transfer 
denotes energy given by the shell Q, while negative transfer corresponds to energy 
received by this shell. In this case, kinetic energy in the shell Q = 20 is mostly re- 
ceived from the shell K = 18 (negative peak), and given to K = 22 (positive peak). 
This function represents the local and direct transfer of kinetic energy to small 
scales. There are no noticeable differences in this transfer between the Hall-MHD 
runs (e = 0.05 and 0.1), and the MHD run (e = 0). 

The curve for early times (kinematic regime) corresponds to the initial exponen- 
tial growth of magnetic energy, and is a time average properly normalized. As time 
evolves and the magnetic energy grows, the amount of kinetic energy transfered 
to small scales diminishes, since a larger amount of kinetic energy at large scales 
is turned into magnetic energy. This effect was previously observed in MHD runs 
(Mininni et al., 2005a). 

Figure 5 shows r|| HD (A', Q = 20), the transfer of magnetic energy at the shell 
Q = 20 to magnetic energy in shells K due to the advection by the bulk velocity 
field. As in the case of Tjju, the transfer is local and the shell Q = 20 receives 
most of the energy from K = 18 (negative peak) and gives energy to the shell 
K = 22 (positive peak). Again, no significant differences are observed between the 
three runs with different values of e, except that this transfer, in amplitude, gets 
substantially stronger as e (and B) increases. 

The total shell-to-shell transfer of magnetic energy is given by T~]^p^ plus . 
Figure 6 shows the T]^g (K, Q) transfer at Q = 20. As in the previous cases, 
positive transfer denotes energy is given from the shell Q = 20 to shells K, while 
negative transfer indicates the shell Q receives energy from K. The "?s|^ u transfer 
is small during the kinematic regime, but grows as the small scales reach nonlinear 
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t igure 6. T^£ U (K, Q = 20) in the kinematic regime (solid line), at t = 26 (dotted line), 
and at t = 45 (dashed line) in the run with e = 0.1. 
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Figure 7. Tbu(K = 20, Q) in the kinematic regime (solid line), at t — 26 (dotted line), 
and at t — 45 (dashed line) in the run with e = 0.1. The inset shows a blow up of the last 
transfer. 

saturation. Although this transfer is noisier than the previous terms studied, two 
regions can be identified at late times. Around Q = K = 20, the transfer is local 
and direct: positive and negative peaks can be observed at K s=s 25 and K « 15, 
indicating energy is received and given respectively by the shell Q from and to 
these wavcnumbers. On the other hand, at large scales (up to K « 10) a region with 
positive transfer can also be identified. This region indicates a non-local and inverse 
transfer of energy: the shells with K between 1 and 10 receive magnetic energy 
from the shell Q = 20. This combination of a local direct transfer of energy and 
a non-local inverse transfer is characteristic of the Hall term, and is in qualitative 
agreement with the turbulent dissipation derived in Sec. 3 where it was shown that 
it can take negative values. 
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Figure 8. Tbu{K,Q — 20) at three different times: the kinematic regime (solid line), 
t — 26 (dotted line), and t = 45 (dashed line) for the run with e = 0.1. 

The remaining transfer term is Tbij{K, Q), which when positive represents trans- 
fer of kinetic energy from the shell Q to magnetic energy in the shell K. Although 
the expression of this transfer function is equal for MHD and Hall-MHD, the trans- 
fer is modified by the Hall currents. The reason for this can be explained in two 
ways. On the one hand, the expression of the a-effect in Sec. 3 is modified by the 
Hall term, and this term represents transfer of energy from the turbulent velocity 
field to the mean magnetic field. On the other hand, waves are expected to give 
non-local coupling between the velocity and magnetic fields (see e.g. Iroshnikov, 
1963; Kraichnan, 1965, in MHD). In Hall-MHD, the non-dispersive Alfven waves 
of MHD are replaced by dispersive circularly polarized waves and as a result, the 
coupling between the two fields should also be modified. 

Figure 7 shows Tbu(K = 20, Q), the energy transfered to the magnetic field at 
K = 20 from the velocity field at shells Q. In the kinematic regime this transfer is 
non-local and similar to the MHD transfer (Alexakis et al., 2005b; Mininni et al., 
2005a): the magnetic field at K = 20 receives energy from the large scale flow at 
Q = 3 and from all turbulent scales up to Q « 20. However, at late times the 
transfer is strongly modified. The magnetic field at K = 20 still receives energy 
from a broad range of wavenumbers Q smaller than K (as was found in Alexakis 
et al., 2005b), but it also receives energy from larger wavenumbers (Q ss 22), and 
gives energy to the velocity field at slightly smaller wavenumbers (Q sa 18). Note 
that this indicates that in Hall-MHD a magnetic field at a given scale can give rise 
to velocity fluctuations at larger scales, a process studied by Mahajan et al. (2005b) 
and referred there as the reverse dynamo. 

Figure 8 shows Tbu(K, Q = 20), the energy received by the magnetic field at all 
wavenumbers K from the velocity field in the shell Q = 20. During the kinematic 
regime, the velocity field in this shell gives energy to all magnetic shells, although 
the transfer peaks at wavenumbers larger than Q. But in the saturated regime, the 
transfer changes drastically again. The magnetic field at wavenumbers K smaller 
than Q « 16, and in shells between 20 and 23 gives energy to the velocity field 
(negative transfer), while the magnetic field in shells between K m 16 to 20 and 
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Figure 9. T^g (K,Q) (left column), and Tbu(K,Q) (right column) at three different 
times: the kinematic regime (top), t — 26 (middle), and t — 45 (bottom) for the run with 
e = 0.1. In all figures, K is on the x axis and Q in the y axis. Shading go from dark 
(T < 0) to light (T > 0). 



for K > 23 receives energy from the velocity field (positive transfer). This is just 
the counterpart of Tbu(K,Q) for constant K, and again shows that in Hall-MHD 
a small scale magnetic field can create large scale flows. 

Figure 9 shows shaded plots of T^g^(K,Q) and Tbu{K,Q) at different times. 
These are the two transfers that are strongly modified by the Hall currents, and 
the figures allow for a study of the terms for all values of K and Q. Although noisy, 
a characteristic pattern can be recognized in Tg^^. ^ s time evolves and the mag- 
netic energy grows, the relative importance of this term grows. For wavenumbcrs 
K-i Q ~ k-Haii ~ 10, the function is positive (light) near and below the diagonal 
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Figure 10. Energy fluxes in the run with e = 0.1 at t = 45: Il£f u (fc) (solid line), 

n| I J ID (fc) (dash-dotted line), Uuu(k) (dashed), and Yl B u(k) (dotted). The thick line is 
the total flux in the simulation. 

K = Q, and negative (dark) near and above this diagonal. This region close to 
the diagonal represents local and direct transfer of energy: a cut at constant Q 
shows that close to the diagonal the shell Q receives energy from neighboring shells 
with K < Q (negative Tj^f^) and gives energy to neighbor shells with K > Q 
(positive 7g^^). As we move far from this diagonal, the sign of the regions above 
and below the diagonal changes. This indicates a non-local and inverse transfer of 
magnetic energy, from small to large scales, in agreement with the expression for 
the turbulent magnetic diffusivity obtained in Sec. 3. 

The Tbu(K, Q) also shows an interesting behavior as a function of time. During 
the kinematic regime, Tbu is positive in a triangle defined by K > Q. This indicates 
that the velocity field in a given shell amplifies the magnetic field in that shell and 
all the shells with larger wavenumber (smaller scales). Also a strong band around 
Q = 3 is observed, indicating that the velocity field in the energy injection band 
gives a lot of energy to the magnetic field. These results are similar to the kinematic 
MHD dynamo (see Mininni et al., 2005a). However, at late times an inverse process 
can be identified close to the diagonal K = Q, Above it, Tbu is positive, while 
below it, it is negative. This represents transfer of magnetic energy from a shell K 
to kinetic energy in slightly smaller wavenumbers Q. 

Since Tj^l 1 ^ and Tbu gi ve both direct and inverse transfers of energy (locally or 
non- locally), it is of interest to quantify which direction wins when all the contri- 
butions to the transfer are added. To this end, we computed the contribution of 
each transfer term to the energy flux. The total energy flux at a wavenumber k is 
given by 



where T = luu + Tb b + Tub + Tbu is the total energy transfer. We can split this 



k 




(5.1) 



K=0 Q 
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flux into the energy flux due solely to the transfer of kinetic energy 

k 

^uu(k) = E Y,Tuu{K,Q), (5.2) 

K=0 Q 

the flux due to the transfer of magnetic energy Hbb = n^^^ + n^J; , where 

nMl iD w = EE^I HD (A',g) J (5.3) 

K=0 Q 

nf^ = EE r *B a V,Q), (5.4) 

K=0 Q 

and the hybrid flux due to interactions between the velocity and magnetic fields 

k 

n BU (k)= Y,Y,[Tbu(K,Q) +T UB (K,Q)] . (5.5) 

K=0 Q 

To compute the fluxes, the transfer functions arc not normalized. 

Figure 10 shows the partial energy fluxes at t = 45. Since all transfer functions 
were computed up to K, Q = 40, the partial fluxes go to zero artificially at this 
wavenumbcr, although in the simulation the total energy flux goes to zero only at 
the maximum resolved wavenumber k max . 

The total flux is positive at wavenumbers larger than ko (the energy injection 
band), indicating a direct cascade of the total energy. At wavenumbers smaller than 
ko the total flux is negative, an evidence of large scale dynamo action. A substan- 
tial portion of the total flux is due to the transfer of energy from the kinetic to 
the magnetic reservoirs (Hsu), an d this contribution to the flux is positive at all 
wavenumbers (larger than the forced ones) indicating a net direct transfer of the 
energy. We note that this flux is due to the non-local Tub and Tbu transfer terms. 
The flux due to the transfer of kinetic energy Tluu is also positive at all wavenum- 
bers. But the flux due to the transfer of magnetic energy Ii% B 11 ^ s oru y positive at 
wavenumbers larger than kuaii- For wavenumbers smaller than k « 10 ~ kjjaii, 
Ilgf 11 chan ges sign, giving as a result a net inverse transfer of magnetic energy, 
from small to large scales. This indicates that a magnetically dominated Hall-MHD 
system could display backscatter of the magnetic energy Magnetic fluctuations at 
small scales could give rise to large scale magnetic fields, as is also implied by the 
expression of (3 found in Sec. 3. 

Note that although in MHD the inverse cascade of magnetic helicity can give a 
similar result, the backscatter predicted in Hall-MHD by the turbulent diffusivity 
is novel, since it can take place even in the absence of helicity in the fields. To 
illustrate this we show results of non-helical magnetically dominated simulations in 
Sec. 6. It is worth noting that at wavenumbers smaller than kuaii, the Hall term 
increases the flux of magnetic energy to smaller scales, thus also in agreement with 
results showing the Hall currents increase the amount of small scale perturbations 
(Birn ct al., 2001; Laveder et al., 2002a,b; Morales et al., 2005). 

5.2. Dependence with e 

Now wc discuss in detail the dependence of the results as e (or the Hall scale) is 
varied. To this end, wc consider the runs with e = 0.1, 0.05, and (MHD). As 
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Figure 11. T BU (K = 20, Q) at t = 45 for e = (solid line), 0.05 (dotted line), and 0.1 

(dashed line). 

previously mentioned, the transfer terms Tjju(K, Q) and Tbb(K, Q) do not show a 
dependence with the amplitude of the Hall effect. These terms give direct and local 
transfer of energy to small scales, as in MHD (Alexakis et al., 2005b). As a result, 
we will discuss the change in the remaining transfer terms as e is varied. 

Since the Hall effect is more relevant when the magnetic field is stronger, we 
will consider the transfer terms at t = 45, when a large scale magnetic field is 
present and the small scales have reached saturation. Since we will examine runs 
with different values of e at the same time, in this subsection the transfers are not 
normalized using the energies. 

The transfer of magnetic energy due to the Hall term T^g-{K,Q = 20) is of 
course zero in the MHD case (e — 0). As e is increased, except for an increase in 
its amplitude (not shown), no significant differences are observed and its behavior 
is similar to the one examined in Sec. 5. 

Figure 11 shows the behavior of Tbu(K = 20, Q) (the transfer of kinetic energy 
in the shells Q to magnetic energy in the shell K = 20) as e is varied. The strong 
peak at Q = 3 is associated with the injection band. This transfer is non-local in 
the three runs, as is evidenced by the positive plateau from Q rs 3 to Q ps 16. As 
a result, the velocity field in all these shells gives energy to the magnetic field at 
K = 20. As e is increased, a local transfer grows in the neighborhood of Q = 20. The 
velocity field at wavenumbcrs K slightly larger give energy to the magnetic field at 
K = 20, while the magnetic field gives energy to the velocity field at wavenumbers 
slightly smaller (Q ps 18). 

We can compute the energy flux as e increases. Since the transfer of kinetic 
energy is not changed, we will focus on two contributions to the total flux: the flux 
of magnetic energy U.BB(k), and the hybrid flux Usuik) due to the terms turning 
kinetic into magnetic energy and vice- versa. The magnetic energy flux Hbb = 
^bb' + I^bb is shown in Fig. 12. At scales larger than kHalU negative flux of 
magnetic energy is observed, giving backscatter of magnetic energy to large scales. 
As e is increased, the amplitude of the backscatter grows, and the wavenumbcr 
where the flux changes sign moves to larger k. 
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Figure 12. n SB (fc) at i = 45 for e = (solid line), 0.05 (dotted line), and 0.1 (dashed 

line). 
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Figure 13. IL B u(k) at t = 45 for e = (solid line), 0.05 (dotted line), and 0.1 (dashed 

line) . 

Figure 13 shows the flux Hgij(k). At wavenumbers smaller than the forcing 
wavenumber (k = 3), the flux is negative. This is a signature of large scale dynamo 
action: the magnetic field at large scales is fed by the small scale velocity field. 
Remarkably, as e is increased, the amplitude of the negative flux at large scales 
increases. This is in good agreement with dynamo simulations where the large scale 
magnetic field was observed to grow faster in the presence of Hall currents (Mininni 
et al, 2003a, 2005b). 

5. 3. Transfer of magnetic helicity 

We discuss briefly the transfer of magnetic helicity in the saturated case t = 45. To 
study the transfer associated with the inverse cascade of magnetic helicity at scales 
larger than the forcing scale, a large separation between this scale and the largest 
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Figure 14. Th(K, Q = 20) normalized by the magnetic helicity at the shell Q = 20, at 
t = 45 for e = (solid line), 0.05 (dotted line), and 0.1 (dashed line). 

scale in the box is needed. At a fixed spatial resolution, this reduces the Reynolds 
numbers, and as a result reduces also the separation between the Ohmic scale and 
the Hall scale. This study is beyond the aim of this work. But we want to point out 
a remarkable feature observed in the transfer of magnetic helicity at scales smaller 
than the forcing scale. 

Figure 14 shows the transfer Th (K, Q = 20) normalized by the magnetic helicity 
in the shell Q = 20 at t = 45. The transfer is mostly local in the three simulations, 
peaking at wavenumbers K slightly smaller and larger than Q. However, as e is 
increased the transfer rate of magnetic helicity is strongly quenched. This slow down 
in the transfer in Hall-MHD explains the behavior observed in Fig. 3. In MHD and 
Hall-MHD dynamos, the external mechanical forcing generates equal amounts of 
magnetic helicity of opposite sign at scales smaller and larger than the forcing band 
(Seehafer, 1996; Brandenburg, 2001; Mininni et al., 2003a). Since the transfer of 
magnetic helicity between different shells in the Hall-MHD runs is almost stopped, 
it takes more time for the magnetic helicity at scales smaller than the forcing scale 
to reach the dissipative scale where it can be destroyed. As a result, both signs of 
magnetic helicity piles up close to the forcing band, decreasing the growth rate of 
net magnetic helicity at scales larger than the forcing scale, and allowing also for 
the possibility for a sign change of the net magnetic helicity. 



6. Backscatter of magnetic energy in Hall-MHD 

The mechanically forced runs discussed in the previous section show negative flux 
of magnetic energy at large scales due to the Hall effect, in agreement with negative 
values of the turbulent diffusivity. This indicates that in a magnetically dominated 
simulation, backscatter of magnetic energy could be observed if the Hall term is 
strong enough. Note that here we are using the word backscatter to refer to this 
transfer of magnetic energy from the small to the large scales. This is done in 
opposition to the usual terminology of inverse cascades, since we have been unable 
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Figure 15. (a) Kinetic energy spectrum Eu(k) and (b) magnetic energy spectrum 
Es(k) at t — 3, for runs with e = (solid line), 0.05 (dotted line), and 0.1 (dashed line). 



to identify any ideal invariant of the Hall-MHD equations cascading inversely with 
constant flux to the large scales. 

To study this scenario, we did three simulations with e = 0, 0.2, and 0.5. The 
kinematic viscosity and magnetic diffusivity were v = r\ = 3.5 x 10^ 2 , and the 
spatial resolution was N 3 = 128 3 . The initial condition was U = B = 0. The 
system was forced with a non-helical and random electromotive force given by a 
superposition of harmonic modes at wavenumbers k = 9 and 10. The phases of the 
force were changed with a correlation time of r = 1.25 x 10~ 2 , and the time step 
was set to At = 1.5 x 10~ 3 . Note that in the absence of magnetic helicity, no inverse 
cascade is expected in the MHD case (see however Lanotte et al., 1999, for cases 
where a large scale shear is present). 

The system was run until reaching a turbulent steady state. All the quadratic 
invariants (with the exception of the total energy) were verified to be small: the 
magnetic helicity fluctuates in the three runs around zero, both the global quantity 
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Figure 16. (a) Kinetic energy spectrum and (b) magnetic energy spectrum at t — 75 for 
runs with e = (solid line), 0.05 (dotted line), and 0.1 (dashed line). 



as well as its spectral density at each individual Fourier shell. Figure 15 shows 
the kinetic and magnetic energy spectrum at early times (t = 3). The shell of 
wavenumbers associated with the external magnetic force is easily recognized in 
the peak in Fig. 15(b). 

As time evolves, an increase in the magnetic energy at wavenumbers smaller 
than the forcing wavenumber is observed. Figure 16 shows the kinetic and magnetic 
energy spectrum at t = 75, when the system has reached a steady state. The three 
runs are dominated by the magnetic energy (note that the peak in the magnetic 
energy spectrum around the forcing band gives the largest contribution to the 
energy). The spectrum of kinetic energy is similar for the three runs, and large 
scale perturbations are observed because of the injection of kinetic energy by the 
Lorentz force. However, the magnetic energy spectrum is strongly modified as e is 
increased. While the spectra of the three simulations peak in the energy injection 
band, the magnetic energy at K = 1 in the run with e = 0.5 is three orders of 
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magnitude larger than in the MHD run. The magnetic energy in all wavenumbers 
smaller than the forcing wavenumbers increases as e is increased. 

The backscatter of magnetic energy is in good agreement with the negative flux 
YIbb observed in the previous section, and the negative turbulent transport coef- 
ficients derived for the magnetically dominated case. Note that an increase in the 
level of the small scale magnetic fluctuations (for wavenumbers smaller than the 
energy injection wavenumbers) is also observed in Fig. 16(b). 

7. Discussion 

In this work we presented energy transfer in Hall-MHD turbulence as obtained 
from numerical simulations. The properties of the spectral transfer is one of the 
building blocks of turbulence theories, and to the best of our knowledge no attempt 
to study transfer and cascades of ideal invariants in this system of equations had 
been attempted before. 

Before proceeding with the discussion of our results, we have to warn the reader 
about a clear limitation of the numerical results presented. As previously mentioned, 
an astrophysics-like scale separation between the box size, the energy injection scale, 
the Hall scale, and the Ohmic dissipation scale is well beyond today computing 
resources. We tested the dependence of our results as e was varied, but no attempt 
was made to change the Reynolds numbers in our simulations. This being said, we 
believe that even under this limitation, an understanding of the transfer of energy 
between different scales is of uttermost importance for the development of a theory 
of turbulence for Hall-MHD or other extensions of magnetohydrodynamic to take 
into account kinetic plasma effects. 

Direct evidence of nonlocality of the energy transfer was observed. While the total 
energy displays a direct cascade to small scales, in the individual transfer terms, 
both directions (toward small and large scales) were identified. Coupling between 
the magnetic and velocity fields is strongly modified by the Hall effect, and a local 
backscatter of energy from the magnetic field to the velocity field at slightly larger 
scales was observed. This behavior can be expected since the Hall term changes 
the nature of the nondispersive MHD Alfven waves, into dispersive and circularly 
polarized waves. As a result, the nonlinear coupling between the two fields is also 
changed. 

Also a nonlocal backscatter of magnetic energy was observed at scales larger 
than the Hall scale. This backscatter was verified in non-helical magnetically forced 
simulations, where the amplitude of the magnetic field at scales larger than the 
forcing scale was observed to grow in the Hall-MHD simulations, but not in the 
MHD run. In some sense, the magnetic field in Hall-MHD being frozen in the ideal 
case to the electron velocity field, couples non-locally both small scales (the current) 
and large scales (the bulk velocity field). 

All these processes can be partially explained considering transport turbulent 
coefficients estimated from MFT. Unlike MHD, the turbulent diffusivity in Hall- 
MHD is not positive definite. In particular, its expression shows that ion-cyclotron 
waves are more likely to produce large values of negative (backscatter) or positive 
(reconnection) turbulent diffusivity than the whistler mode. 

The transfer of magnetic helicity at small scales was also observed to be quenched 
by the Hall effect. While the mechanisms generating magnetic helicity in the Hall- 
MHD dynamo are the same as in MHD (Mininni ct al., 2003a), the transport of 
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helicity is expected to be changed by the Hall currents (Ji, 1999). As a result of 
the slow down in the transfer rate of magnetic helicity by the Hall effect, the late 
time evolution of the system is not characterized by a maximally helical large scale 
magnetic field as in the MHD case (Pouquet et al., 1976; Mcncguzzi et al., 1981; 
Brandenburg, 2001). 

The Hall term gives a direct transfer of magnetic energy at scales smaller than the 
Hall scale, and an inverse transfer at scales larger than the Hall scale. This finding 
sheds light into the conflicting results reported in the literature, where the Hall 
effect was observed to increase the amount of small scales and magnetic dissipation 
in some cases, and to help large scale reorganization processes in other cases, as 
mentioned in the introduction. 

As a result of this dual direction of the Hall transfer, a change in the power law 
followed by the total energy spectrum can be expected close to the Hall wavenum- 
ber. Steepening of the energy spectrum for wavenumbers smaller than kuaii was 
observed in 2.5D simulations with strong magnetic fields imposed, when the cross- 
correlation between the velocity and magnetic fields was significant (Ghosh et al., 
1996). In three dimensional dynamo simulations where the cross correlation is in 
general small, no change was observed (Mininni et al., 2005b), although a faster 
growth of the large scale magnetic field was found. Given the nonlocal nature of the 
transfer in Hall-MHD, and the scale separation needed to observe a clear change in 
the energy spectrum, probably a huge increase in the spatial resolution is needed 
to confirm it. 
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